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Abstract 

This  paper  provides  a  detailed  review  of  renewably  driven  hydrogen  systems  and  modeling 
approaches  applicable  to  these  systems  that  have  been  reported  over  the  last  two  decades.  Several 
renewable  energy  technologies,  including  solar  photovoltaic,  wind,  and  hydro,  have  been  considered 
as  the  power  source  in  these  reports.  Storage  is  an  important  aspect  of  hydrogen  systems,  and  options 
for  this  are  summarized  here.  Utilization  systems  include  fuel  cells  as  well  as  a  variety  of  thermal  uses. 
Some  of  the  reported  studies  have  addressed  residential  applications  whereas  others  were  related  to 
commercial  scale  systems.  This  paper  particularly  emphasizes  aspects  of  modeling  of  the  various 
components  for  the  renewable  hydrogen  system.  Based  on  the  literature  on  this  area,  conclusions  are 
provided  on  the  current  understanding  as  well  as  future  work  related  to  this  topic. 
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1.  Introduction 

As  concerns  grow  about  degradation  of  the  environment,  availability  of  conventional 
fuels,  and  the  political  viability  of  nuclear  power,  alternatives  for  energy  supplies 
need  to  be  examined.  One  of  the  hopes  for  the  future  is  that  renewably  based  approaches 
for  energy  can  be  used  as  replacements  for  our  currently  more-widely  used  sources. 
However,  many  of  the  renewable  energy  sources  have  some  degree  of  temporal 
output,  causing  problems  for  power-on-demand  requirements.  For  this  reason, 
some  form  of  storage  will  undoubtedly  be  of  value.  Also,  there  will  be  some  future 
requirements  for  a  transportable  fuel  that  may  or  may  not  be  able  to  be  satisfied 
by  batteries. 

One  approach  to  be  considered  is  use  of  hydrogen  as  a  significant  element  in  our  energy 
future.  Major  advantages  of  using  hydrogen  include:  firstly,  it  is  a  clean  source  of  energy; 
secondly,  it  can  be  used  for  seasonal  storage  of  energy;  thirdly,  it  can  be  transported  (say  as 
an  aircraft  fuel),  and  finally  it  can  be  used  to  produce  electrical  as  well  as  thermal  forms  of 
energy. 

The  concept  of  integrating  renewable  energy  with  hydrogen  systems  was  given  some 
serious  consideration  in  the  1970s  [1,2].  Fischer  [3]  provided  an  interesting  review  on 
experimental  and  analytical  studies  of  early  photovoltaic/electrolyzer  systems  used  for 
hydrogen  production.  Most  of  the  experimental  studies  reported  were  on  a  laboratory  scale. 

In  last  three  decades,  intense  research  has  been  conducted  around  the  world  on  this 
topic.  Emphases  are  not  only  to  improve  the  performance  of  existing  hydrogen  production, 
storage  and  utilization  technologies,  but  also  to  integrate  them  effectively  with  renewable 
energy  sources.  As  noted  above,  one  of  the  key  elements  of  this  is  that  it  provides  a  good 
buffer  to  time- varying  energy  inputs. 

Fig.  1  shows  a  generic  renewable  hydrogen  system.  The  renewable  energy  source  can  be 
of  any  type,  including  time-varying  forms  such  as  solar,  wind,  and  hydroelectric.  Hydrogen 
can  be  produced  by  electricity  from  these  kinds  of  systems  as  well  as  by  chemical  and 
biological  means  [4]. 

Storage  is  an  important  aspect  of  this  type  of  system.  This  is  particularly  true  if  this  is 
driven  by  a  time-varying  energy  source  but  it  is  also  needed  if  there  is  a  widely  varying  time 
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Nomenclature 


Apv  area  of  PV  module  (m  ) 

a  modified  ideality  factor  (< a  =  npvkTpvNs/q ) 

ac ,  bc  empirical  constants  for  cell  temperature  calculation 

Afc  active  area  of  cell  (cm  ) 

age  operation  time  for  fuel  cell  (h) 

aref  reference  modified  ideality  factor 

c*H+  proton  concentration  at  the  cathode  membrane/gas  interface  (mol/cm3) 

<fHi  liquid-phase  concentration  of  hydrogen  at  the  anode/gas  interface 
(mol/cm3) 

ch,o  water  concentration  at  the  cathode  membrane/gas  interface  (mol/cm  ) 

Cq ,  concentration  of  oxygen  at  cathode  membrane/gas  interface  (mol/cm3) 

Ccw  thermal  capacity  of  cooling  water  (J/K) 

Cmh  effective  heat  capacity  (J/°C) 

Cp  specific  heat  at  constant  pressure  (J /mol  K) 

Ct,el  overall  heat  capacity  of  electrolyzer  (J/K) 

F  Faraday’s  number  (96485) 

Gnoct  instantaneous  solar  irradiation  at  SRC  (W/m  ) 

Gt  instantaneous  total  solar  irradiation  on  PV  (W/m  ) 
i  current  density  draw  from  fuel  cell  (A/cm  ) 

I ei  electrolyzer  current  (A) 

ll,  ref  reference  light  current  (A) 

IL  light  current  (A) 

iL  maximum  current  density  (A/cm  ) 

lm, ref  reference  maximum  power  point  current  of  PV  cell  (A) 

Im  maximum  power  point  current  of  PV  cell  (A) 

Io,ref  reference  diode  reverse  saturation  current  (A) 

I0  diode  reverse  saturation  current  (A) 

Ipv  current  of  PV  cells  (A) 

I  sc, ref  reference  short  circuit  current  (A) 

I sc  short  circuit  current  (A) 

k  Boltzmann’s  constant  (1.381  x  10-23  J/K) 

k°cr  k°c  reaction  rate  constants  for  anode  and  cathode,  respectively  (cm/s) 
kceii  empirical  term  to  represent  the  apparent  rate  constant  for  anode  and 
cathode 

kDR  empirical  term  to  represent  the  change  in  catalytic  activity  with  age  (V/hK) 
keff  effective  heat  loss  coefficient  (W/°C) 

L  membrane  thickness  (cm) 

LMTD  log  mean  temperature  difference  for  heat  exchanger  (K) 

n  number  of  electrons  transferred  per  reaction  (2) 

nc  number  of  electrons  involved  in  the  reaction 

Nc>ei  number  of  cells  in  electrolyzer  stack 

Ncfc  number  of  cells  in  fuel  cell  stack 

ngas  hydrogen  flow  rate  (mol/s) 

vlh 2  molar  flow  rate  of  hydrogen  (mol/s) 
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fiH2o  molar  flow  rate  of  water  (mol/s) 

ho '2  molar  flow  rate  of  oxygen  (mol/s) 

ripoiy  poly  tropic  index 

npv  parameter 

Ns  number  of  cells  in  series 

nst  number  of  moles  of  hydrogen  inside  the  storage  tank  (mol) 

p  operating  pressure  of  electrolyzer  cell  (N/m  ) 

Pa,  ht  actual  available  power  of  hydro  turbine  (W) 

Pa,  wt  actual  available  power  of  wind  turbine  (W) 

<y 

pc  in  pressure  at  the  outlet  of  compressor  (N/m  ) 

<y 

Pc, out  pressure  at  the  outlet  of  compressor  (N/m  ) 

P comp  power  required  by  compressor  (W) 

peq  equilibrium  pressure  of  hydrogen  gas  (N/m2) 

P h2  partial  pressure  of  hydrogen  (N/m  ) 

P 02  partial  pressure  of  oxygen  (N/m  ) 

*y 

Pref  reference  pressure  (N/m  ) 

pst  pressure  inside  the  storage  tank  (N/m  ) 

Pth,  wt  theoretical  available  power  of  wind  turbine  (W) 
q  electronic  charge  (1.609  x  1(T19  As) 

Qcooi  heat  lost  to  cooling  media  (W) 

Qgen  heat  generated  by  elecrolyzer  (W) 

Qioss  heat  loss  (convection  and  radiation  to  ambient)  by  electrolyzer  (W) 
Qmh  total  heat  production/consumption  (W) 

R  gas  constant 

r i,  r2,  s,  1 1,  t2,  h  empirical  constants  for  electrolyzer 
R electronic  resistance  to  qow  0f  electrons  at  the  electrodes  (Q) 

rm  membrane  specific  resistivity  for  the  flow  of  hydrated  protons  (Q  cm) 

Rproton  resistance  to  the  flow  of  protons  in  the  cell  (Q) 

Rs  ref  reference  series  resistance  (Q) 

Rs  series  resistance  (Q) 

Rsh,ref  reference  shunt  resistance  (Q) 

Rsh  shunt  resistance  (Q) 

Rt  ei  overall  thermal  resistance  of  electrolyzer  (W/K) 

Spv>ref  reference  effective  absorbed  solar  irradiation  (W/m2) 

Spv  effective  absorbed  solar  irradiation  (W/m  ) 

T  operating  cell  temperature  (K) 

T a, n oct  ambient  temperature  at  SRC  (K) 

Ta  ambient  temperature  (K) 

Tcw>i  cooling  water  temperature  at  inlet  (K) 

Tcw>0  cooling  water  temperature  at  outlet  (K) 

Tm,pv  back  surface  temperature  of  PV  module  (K) 

Tnoct  cell  temperature  at  standard  rated  condition  (SRC)  (K) 

Tpv  cell  temperature  of  PV  (K) 

Pref  reference  temperature  (K) 

C0  standard  cell  potential  (V) 

U0  fc  standard  cell  potential  for  fuel  cell  (V) 
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UAhx  product  of  overall  heat  transfer  coefficient  and  area  of  for  heat  exchanger 
(W/K) 

Uei  cell  potential  for  electrolyzer  (V) 

Ufc  overall  cell  potential  for  fuel  cell  (V) 

Ul,noct  overall  heat  transfer  (convection  and  radiation)  coefficient  of  PV  at  SRC 
(W/m2  K) 

UL  overall  heat  transfer  (convection  and  radiation)  coefficient  of  PV  (W/m  K) 
Um>ref  reference  maximum  power  point  voltage  of  PV  cell  (V) 

Um  maximum  power  point  voltage  of  PV  cell  (V) 

Uoc,ref  reference  open-circuit  voltage  (V) 

Uoc  open-circuit  voltage  (V) 

Upv  voltage  of  PV  cells  (V) 

Urev  reversible  cell  potential  (V) 

Utn  thermo-neutral  potential  (V) 

Vst  volumetric  capacity  of  storage  tank  (m  ) 

WS  wind  speed  (m/s) 

a,  ft,  y,  5,  £  empirical  constants  for  calculation  of  specific  heat  at  constant  pressure 
otc  chemical  activity  parameter  for  the  cathode 

AG  change  in  Gibbs  free  energy  (J/mol) 

A  Ge  standard  state  free  energy  of  the  cathode  reaction  (J/mol) 

A  Gec  standard  state  free  energy  of  chemisorption  from  the  gas  state  (J/mol) 

AH  change  in  enthalpy  of  formation  (J/mol) 

A H°  standard  (105N/m2  and  298  K)  enthalpy  of  formation  (J/mol) 

AS  change  in  entropy  of  formation  (J/mol  K) 

A S°  standard  (105N/m2  and  298  K)  entropy  of  formation  (J/mol  K) 

A  Tpv  temperature  difference  between  cell  and  module  back  surface  at  a  1000 

W/m2  irradiance  (°C) 

Oact  activation  losses  (V) 

f]c,pv  module  efficiency  to  convert  incident  solar  irradiation  to  electrical  energy 
rjcomp  compressor  efficiency 

Oconc  concentration  losses  (V) 

rjohm  ohmic  losses  (V) 

£i,  (^2?  £3,  ^4  empirical  parameters  for  fuel  cell 


variation  of  use.  This  can  be  accomplished  using  compressed  gas  storage  tanks,  liquefied 
hydrogen  storage,  metal  hydrides  or  carbon-based  materials.  We  also  consider  short-term 
electrical  storage  as  an  option.  This  could  be  batteries  or  grid  connection,  if  the  latter  is  available. 


Fig.  1.  Generic  renewable  hydrogen  system. 
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Loads  for  houses  can  be  classified  on  an  electrical  or  thermal  basis.  Electrical  loads  can 
be  satisfied  by  fuel  cells  (stationary  and  automobile  application)  whereas  hydrogen  can  be 
combusted  for  satisfying  thermal  loads,  e.g.,  IC  engines  (cars),  cooking,  water  heating, 
HVAC  systems,  etc. 

A  variety  of  studies  that  have  been  performed  on  combined  renewable  energy  and 
hydrogen  systems  can  broadly  be  classified  as  follows. 

1.1.  Renewable  hydrogen  generation  system  (RHGS) 

The  system  generates  hydrogen  from  renewable  energy  sources  and  can  be  either 
stand  alone  or  grid  connected.  This  system  is  a  typical  hydrogen  refueling 
station.  Note  that  for  residential  applications  such  a  system  can  fulfill  the  hydrogen 
requirements  for  fuel  purposes  only.  A  description  has  been  given  for  the  initial 
design,  operation  and  cost  optimization  of  the  system  with  photovoltaic  cells  connected 
to  an  electrolyzer  [5-9].  Experimental  data  on  small-scale  systems  with  proton 
exchange  membrane  (PEM)  electrolyzers  have  been  summarized  [10].  Experimental 
data  for  RHGS  with  PV  that  was  installed  on  a  house  in  Switzerland  have  been 
reported  [11]. 

One  of  the  first  studies  including  description  of  component  models  to  simulate  the 
performance  of  the  RHGS  was  given  by  Carpetis  [12].  Others  have  described  simulations 
performed  on  the  TRNSYS  software  for  similar  systems  [13].  Hydrogen  production  as  a 
function  of  maximum  power  point  and  tilt  angle  has  been  given.  The  performance 
modeling  of  the  RHGS  (with  PV)  was  reported  [14].  Further,  cost  optimization  studies  of  a 
photovoltaic  electrolyzer  system  using  HYSOLAR  along  with  performance  simulations 
were  reported  by  Seigel  and  Schott  [15].  Other  cost  optimization  studies  have  also 
appeared  in  the  literature  [16,17]. 

Some  of  the  fairly  recent  studies  discussed  the  preliminary  system  design  and  economic 
feasibility  of  a  combined  PV  and  PEM  electrolyzer  [18,19].  Bilgen  developed  models  to 
study  performance  and  cost  optimization  of  PV  electrolyzer  systems  [20].  The  paper 
reported  comparisons  for  fixed  and  tracking  panels.  Others  have  reported  the  optimum 
hydrogen  cost  for  a  small-scale  (<  100  W)  directly  connected  RHGS  with  PV  [21].  Another 
study  has  focused  on  the  use  of  PSCAD/EMTDC  software  for  the  design  and  simulation 
of  transient  performance  of  RHGS  with  PV  [22]. 

1.2.  Renewable  hydrogen  generation  and  utilization  system  (RHGUS) 

This  system  not  only  generates  hydrogen  from  renewable  energy  sources,  but  also 
includes  a  load  that  utilizes  hydrogen.  The  load  this  system  serves  could  be  electrical  and / 
or  thermal  in  nature.  The  system  can  be  either  stand  alone  or  grid  connected.  Such  a 
system  includes  on-site/off-site  hydrogen  generation  also.  The  latter  refers  to  a  system  in 
which  hydrogen  is  generated  at  a  renewable  source  facility  and  transported  by  some  means 
to  the  point  of  use. 

Design  considerations  based  on  efficiencies  and  energy  calculations  for  various 
components  of  RHGUS  with  PV  as  an  energy  source  were  reported  [23,24].  The  long¬ 
term  operating  experience  for  large-scale  demonstration  plants  of  RHGUS  systems 
installed  in  Germany  has  been  described  [25,26].  Griebhaber  summarized  the  energy 
consumed  and  generated  for  RHGUS  installed  on  self-sufficient  Solar  House  [27].  The 
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efficiencies  for  various  components  of  RHGUS  that  supplied  continuous  power  to  an  air 
compressor  have  been  noted  [28].  Experimental  data  for  the  performance  of  a  combined 
PV  and  wind  turbine-operated  RHGUS  for  remote  applications  were  summarized.  A 
programmable  load  was  used  to  provide  the  energy  demand  for  the  house  [29].  Others 
reported  the  performance  of  RHGUS  with  PV  installed  to  supply  continuous  power 
(200  W)  for  a  remote  application.  A  metal  hydride  bed  was  used  to  store  the  hydrogen 
generated  by  an  electrolyzer  [30]. 

Models  were  developed  and  overall  system  simulations  were  performed  to  optimize  the 
size  of  PV  array  area  for  the  RHGUS  system  [31].  The  program  H2PHOTO  was 
developed  and  simulations  were  performed  based  on  simple  models  (electrical  and 
thermal)  to  design  and  predict  the  performance  of  1.3  kW  RHGUS  system  [32].  TRNSYS 
was  used  to  simulate  the  annual  performance  of  RHGUS  system  for  residential 
purposes.  A  catalytic  burner  was  used  for  providing  heat  to  the  house  but  the  sizing  of 
the  components  depends  on  both  the  electrical  and  thermal  loads  of  the  house  [33]. 
Vosen  and  Keller  [34]  described  time  dependent  annual  simulations  for  the  RHGUS  with 
PV  for  stand-alone  applications  using  state  of  charge  and  neural  network  control 
systems.  The  cost  comparison  between  hydrogen  only,  battery  only  and  combined 
system  was  also  provided.  Others  developed  a  model  to  estimate  hourly  solar  radiation  on 
a  tilted  surface  from  monthly  average  daily  insolation.  This  was  further  used  to 
estimate  the  energy  available  from  a  PV  (model  based  on  efficiencies)  and  wind  turbine 
combined  power  source  for  a  RHGUS  [35].  The  RHGUS  system  for  stand-alone 
residential  application  was  simulated  for  PV  only,  PV-micro-hydro  turbine  and 
combined  PV,  wind  and  micro-hydro  turbine  as  energy  sources  [36,37].  A  thermo- 
economic  analysis  was  performed  for  RHGUS  with  PV  [38].  The  optimization  to 
minimize  the  investment  cost  of  RHGUS  with  PV  and  micro-hydro  turbine  was  carried 
out  [39].  The  results  gave  local  minima  but  it  was  not  sure  that  the  model  gave  global 
minima.  Short-term  performance  of  the  RHGUS  with  PV  and  wind  turbine  for  the 
control  system  and  power  conditioning  devices  developed  by  authors  was  also  reported 
[40].  For  the  similar  system,  the  effect  of  operation  of  electrolyzer  with  oxygen  recovery 
on  the  performance  of  overall  system  has  been  described  [41].  The  dynamic  models  of 
various  components  of  RHGUS  with  PV  and  wind  turbine  were  developed  and 
compared  with  experimental  results.  Further,  these  models  were  used  to  simulate 
the  system  performance  for  residential  purposes  [42].  The  importance  of  the  control 
strategies,  based  on  simulations  performed  on  TRNSYS,  to  operate  the  RHGUS 
with  PV  has  been  noted  [43].  Annual  simulation  for  RHGUS  with  PV  and  wind 
turbine  was  performed  for  an  average  load  of  1  kW  using  Hybrid2  software  developed 
by  National  Renewable  Energy  Faboratory  (NREF)  [44].  Performance  simulation 
of  the  RHGUS  system  with  PV  and  reformer,  installed  in  Stockholm,  Sweden,  to 
obtain  hydrogen  from  biogas  has  been  outlined  in  the  literature  [45].  Shapiro  et  al. 
reported  the  performance  of  an  electrolyzer  (up  to  200  psig)  placed  inside  the  storage  tank 
that  can  be  combined  with  RHGUS  [46].  Dynamic  modeling  of  hydrogen  side 
(electrolyzer,  storage,  and  fuel  cell)  of  RHGUS  has  been  performed  for  application  for 
residential  application  when  combined  with  PV  array  [47].  The  annual  performance 
modeling  of  RHGUS  with  PV  for  a  standard  residential  house  in  Fas  Vegas,  NV  has  been 
performed  by  Deshmukh  and  Boehm  [48].  The  hydrogen  needed  for  fuelling  the 
automobile  and  other  residential  usage  (water  heating,  clothes  drying,  cooking,  etc.)  has 
been  considered. 
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2.  Modeling  by  component 

This  section  reviews  various  approaches  to  modeling  of  various  aspects  of  renewable 
hydrogen  systems.  Models  considered  are  based  upon  their  applicability  for  system 
simulation  purposes.  Both  components  and  input  phenomena  are  included  here. 

2.1.  Input  resources  and  hourly  data 

Modeling  of  various  components  requires  a  variety  of  input  data.  Hourly  averaged 
annual  data  for  solar  irradiation,  temperature,  wind  speed,  wind  direction,  and  water  flow 
are  generally  needed  for  modeling. 

The  hourly  averaged  solar  irradiation  data  for  US  can  be  found  on  Renewable  Resource 
Data  Center  developed  by  NREL  [49].  Wind  data  can  be  downloaded  from  the  NREL 
website  [50].  Monthly  wind  rose  chart  information  averaged  over  thirty  years  can  be 
obtained  from  United  States  Department  of  Agriculture  website  [51].  This  can  be  helpful 
for  the  design  of  wind  turbine  applications  for  a  particular  location.  For  annual 
performance,  hourly  wind  speed  and  direction  data  are  needed  for  a  particular  location. 
Modeling  of  hydro  turbines  needs  hydrograph  data  which  are  listed  on  the  United  States 
Geological  Survey  website  [52].  This  data  can  help  to  design  the  micro-hydro  turbines  and 
predict  their  annual  performance. 

Hourly  data  obtained  can  be  further  used  to  calculate  the  amount  of  solar  irradiation 
incident  on  the  PV  panels.  The  amount  of  incident  irradiation  also  depends  on  the  type  of 
mounting  system  used  for  PV  panels.  Sun-tracking  systems  generally  have  higher  incident 
energy  compared  with  that  on  a  fixed  system  planar  system.  Tracking  systems  can  be 
classified  as  follows: 

•  Single  axis  tracking: 

O  horizontal  east-west, 

O  horizontal  north-south, 

O  vertical. 

•  Two-axis  tracking. 

The  amount  of  energy  available  to  a  variety  of  collector  types  and  orientation  has  been 
developed  for  239  physical  locations  in  the  US  by  the  NREL  [53].  In  doing  this,  hourly 
data  averaged  over  30  years  were  used.  Hourly  data  are  available  for  simulation  purposes. 
One  comprehensive  source  is  the  TMY2  data  [49]. 

Sun-tracking  systems  can  enhance  the  collection  yield  significantly.  When  a  large 
portion  of  the  radiation  is  direct  (compared  with  diffuse),  concentrating  systems  can  be 
used.  The  relation  to  calculate  the  amount  on  incident  solar  irradiation  based  on  type  of 
collector,  whether  tracking  or  not,  is  explained  in  a  number  of  sources  including  the 
excellent  text  by  Duffle  and  Beckman  [54]. 

2.2.  PV  modeling 

Photovoltaic  cells  are  currently  available  in  single  and  multi-junction  forms.  A  problem 
inherent  in  photovoltaic  devices  is  that  they  are  able  to  interact  only  with  a  portion  of  the 
solar  radiation  distribution.  Multi-junction  cells  are  able  to  capture  a  wider  portion  of  the 
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spectrum,  and  they  show  higher  conversion  efficiencies  from  sunlight-to-electricity  than  do 
single  junction  cells.  Commonly  available  photovoltaic  panels  use  single  crystalline, 
polycrystalline,  thin  film  or  multi -junction  cells.  Many  researchers  have  studied  renewable 
hydrogen  systems  with  photovoltaic  cells  (as  mentioned  in  previous  sections)  as  a  source  of 
energy  for  generating  hydrogen. 

Modeling  PV  arrays  for  renewable  hydrogen  systems  can  be  accomplished  by  using  the 
single  diode  model  with  equivalent  circuit  given  in  Fig.  2  [54].  This  model  uses  the  data 
provided  by  manufacturers  (Isc,  Uoc ,  Imp,  Ump,  rIsc,  and  bUoc)  to  model  the  PV 
performance. 

The  IV  curve  is  calculated  as  a  function  of  incident  solar  irradiation  and  cell  temperature 
using  Eq.  (1)  as  follows: 


Io 


exp 


U pi)  T  I uijR, 


pV^S 


a 


C Pv  T  IpvRs 


where  light  current  (/^),  diode  reverse  saturation  current  (I0),  series  resistance  (Rs), 
shunt  resistant  (Rsh),  and  the  parameter  a  are  the  five  characteristics  that  are  needed  to 
solve  Eq.  (1).  Based  on  the  data  from  manufacturers,  Eqs.  (2)— (7)  are  solved 
simultaneously  to  obtain  the  reference  values: 

At  the  short  circuit  condition,  Ipv  =  Isc  and  Upv  =  0, 


sc, ref 


I L, ref  I O, ref 

I  sc, ref  Rs, ref 


exp 


sc,refRs,ref 

CLref 


-  1 


R 


sh,ref 


At  the  open-circuit  condition,  Ipv  =  0  and  Upv  =  Uoc, 


b  —  1 L,ref 


U  oc,ref 

-  1 

exp 

CLref 

U  oc,ref 

Rsh,ref 


Fig.  2.  The  single  diode  equivalent  circuit  model. 
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At  the  maximum  power  point, 

U  m  +  ImRs,ref 


I  m  —  I 


L,ref  I O, ref 


exp 


Clref 


U  it 7  +  I  m  Rsyef 
Rsh,ref 

Also  at  maximum  power  point. 


d  P 


pv 


d  U 


pv 


—  0  —  I  m  +  U  r 


d  / 


pv 


m 


dU 


pv 


m 


-  1 


where 


(4) 

(5) 


d  / 


pv 


d  U 


pv 


1  Q,ref 
U  ref 


exp 


U  m,ref+Im,ref  Rs,ref 
U  ref 


R 


sh,ref 


m 


l  + 


r  /? 

0,ref  t^s/ef 

U  ref 


exp 


E m,ref  El m,ref  Rs,ref 
U  ref 


+ 


Rs,ref 


(6) 


R 


sh,ref 


And,  checking  the  condition  for  estimation  of  the  open-circuit  voltage  predicted  by 
model,  it  is  given  as 


dU0C  U oc(T pv)  -  U oc(T p„>re/) 

=  ~  — 


67 


T  —  T 

±  nil  i 


(7) 


pv  pv,ref 

The  reference  values  are  used  to  calculate  the  five  parameters  at  operating  condition 
using  Eqs.  (8)— (1 3): 


T 


ci  =  aref 


pv 


T 


(8) 


pvref 


II  —  ^  U  Lref  +  Risc  ( T pv  TpVref)\ , 

^  pv,ref 

i  =  i  {  Tpv 

O  O  ref  I  rji 

\A  pVref 

Here, 

Eg  =  Egref[  1  —  C(T pv  —  TpV,ref)\ , 


exp 


'9 


kT 


'9 


kT 


cref 


(9) 

(10) 


.S’ 


pv,ref 

s 


pv 


Finally,  using  the  5  parameters  at  the  operating  conditions,  we  can  either  calculate  the 
current  at  a  given  voltage  or  use  (14)  and  (15)  for  calculating  current  and  voltage  at 
maximum  power  point: 


km  —  It  I 


L 


o 


exp 


Um  T  f  m  Rs 


a 


Um  T  ImRs 


(14) 
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0  —  im  +  u 


m 


ifexp 


U  m~ l~  ImR, 

a 


R 


sh 


1  +  exp 


Bjri+JjnJk 

a 


1  + 


Rx 


R 


sh. 


Performance  of  the  PV  cell  is  very  much  dependent  on  the  cell  temperature  which  itself  is 
a  function  of  solar  irradiation,  ambient  temperature,  wind  speed,  wind  direction  type  of 
mounting,  and  other  factors.  An  energy  balance  for  the  PV  cells  can  be  used  to  find  the  cell 
temperature  from  Eq.  (16)  [54]: 

Tpy  —  T  a  _  Gt  U  l,noct 

T  NOCT  —  T  a,NOCT  G^OCT  Uy 


where 


rj 


c,pv 


InM 


m  ^  m 


ApvGr 


Eq.  (16)  can  be  solved  in  an  iterative  way  to  obtain  the  cell  temperature.  The  cell 
temperature  can  also  be  solved  as  an  explicit  function  of  incident  solar  irradiation,  ambient 
temperature,  and  wind  speed  using  Eqs.  (18)  and  (19)  [55]: 


=  T 


m,pv 


+ 


T 


TNOCT 


where  Tm  is  the  back  surface  temperature  of  module  and  is  given  by 


Ta  +  Gpe 


0 ac+bcWS ) 


Figs.  3-6  show  the  performance  of  PV  panels  of  65  W  rating.  The  data  available  from 
the  manufacturer  are  Isc  =  4.5V,  Uoc  =  2\A,  Im  =  3.95,  Um  16.5,  fiIsc  —  0.00026, 
{Ajjoc  =  —0.85,  Apv  =  0.633,  and  Ns  =  36,  as  reported  in  [54].  Using  the  5-parameter 
model  noted  above,  the  performance  of  the  PV  panels  has  been  calculated.  Fig.  3  indicates 
the  I-U  curves  at  different  levels  of  solar  irradiation  at  a  cell  temperature  40  °C.  The 
maximum  power  point  curve  is  also  shown  in  Fig.  3  .I-U  curve  and  the  power  vs.  U  curves 


'■y 

Fig.  3.  I-U  curve  and  maximum  power  point  curve  for  different  incident  solar  irradiation  (W/m  )  at  40  °C  for  a 
PV  panel. 
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Fig.  4. 


70 

60 

50 

40 

30 

20 

10 

0 


7 

I-  U  curve  and  power  vs.  U  curve  for  a  PV  panel  at  1000  W/m  solar  irradiation  and  25  °C  temperature. 


Voltage  (Volts) 


Fig.  5.  The  effect  of  temperature  on  the  performance  of  a  PV  panel. 


Temperature 


Fig.  6.  Effect  of  temperature  on  maximum  power  of  a  PV  panel. 
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Table  1 


Values  of  empirical  constants  ac,  bc,  and  A Tpv  for  various  types  of  PV  panels  and  mounting  arrangements  [55] 


Module  parameters 

Module 

ac 

bc 

ATpv  (°C) 

Glass/cell/glass 

Open  rack 

—3.47 

-0.0594 

3 

Glass/cell/glass 

Close  roof  mount 

-2.98 

-0.0471 

1 

Glass/cell/polymer  sheet 

Open  rack 

-3.56 

-0.075 

3 

Glass/cell/polymer  sheet 

Insulated  back 

-2.81 

-0.0455 

0 

Polymer/thin  hlm/steel 

Open  rack 

-3.58 

-0.113 

3 

22X  linear  concentrator 

Tracker 

-3.23 

-0.13 
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are  shown  in  Fig.  4  at  a  solar  irradiation  of  1000  W/m  and  cell  temperature  of  25  °C.  The 
effect  of  temperature  on  the  I-U  characteristics  of  the  PV  panel  is  depicted  in  Fig.  5.  The 
effect  of  temperature  on  the  maximum  power  point  for  different  levels  of  solar  irradiation 
is  shown  in  Fig.  6.  It  can  be  observed  that  performance  changes  significantly  as  the  cell 
temperature  changes.  It  is  crucial  to  develop  the  methods  to  keep  the  cell  temperature  as 
low  as  possible  to  improve  the  performance  of  the  PV  cells  (Table  1). 

Other  formulations  for  modeling  photovoltaic  devices  have  been  reported.  For  example, 
King  [55,56]  and  Luft  et  al.  [57]  have  given  examples  of  this.  King  et  al.  [55,56]  developed 
an  empirical  model  for  simulating  the  performance  of  various  PV  technologies  (the  models 
take  electrical,  thermal,  solar  spectral  and  optical  effects  into  consideration). 

Empirical  parameters  found  were  specific  for  each  type  of  module  and  a  list  of  these  is 
maintained  at  Sandia  National  Laboratories  website  [58].  Luft  et  al.  [57]  developed  a 
model  for  predicting  the  current  produced  based  on  the  Isc,  7m,  Um,  Uoc.  Hadj  et  al.  [59] 
developed  different  relations  to  evaluate  the  five  parameters  used  in  Eq.  (1)  for  modeling 
PV  performance.  Comparisons  of  these  models  for  representing  PV  performance  are  given 
in  [60]. 


2.3.  Wind  machine  modeling 


Studies  of  combined  PV  and  wind  power  have  been  reported  [26,29,35,36,40-42,44]. 
Modeling  was  generally  based  on  the  machine  performance  data  available  from  the 
manufacturer.  Fig.  7  shows  the  experimental  and  calculated  (using  (21))  wind  speed  vs. 
power  generated  for  a  2.5  kW  rated  wind  turbine  [61].  With  the  wind  data  for  a  location 
known,  the  amount  of  power  generated  can  be  estimated  over  a  year’s  period. 

Theoretical  available  wind  power  for  a  particular  location  can  be  calculated  using 

Pth,wt  =  \PairA  It  WS 3  •  (20) 

The  power  in  the  wind  must  be  modified  by  the  “fractional  efficiency”  to  determine  the 
power  out  of  the  wind  machine.  Typical  values  of  this  parameter  will  be  in  the  range 
0.15-0.30. 

A  polynomial  fit  can  be  obtained  to  obtain  power  as  a  function  of  wind  speed  for  any 
location  and  a  particular  wind  turbine.  The  result  for  the  data  shown  in  Fig.  7  is  as  given  in 
Eq.  (21): 


Wt  — 


0.177  WS4  +  3.164  WS 3  +  6.757  WS2  -  25.042  1US  +  2.119. 


(21) 
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Fig.  7.  An  example  trace  of  power  output  for  a  range  of  wind  speeds  is  shown  (after  [61]). 


Water  flow  rate  (m3/s) 


Fig.  8.  Maximum  power  generated  for  a  particular  hydro  turbine  is  shown  as  a  function  of  water  flow  rate  and  net 
head  [62], 


2.4.  Micro-hydro  turbine  modeling 

Studies  of  micro-hydro  turbine  powering  hydrogen  generation  have  been 
reported  [36,37,39].  The  modeling  was  done  based  on  the  data  available  from 
the  manufacturers  of  the  equipment  used.  Generally  for  residential  applications, 
low  head  hydro  turbines  (e.g.,  reaction  turbines)  can  be  used.  An  example  of  data 
available  from  a  manufacturer  for  a  reaction  turbine  with  rated  power  2.2  kW  [62]  is  shown 
in  Fig.  8. 

Actual  power  generated  for  the  hydro-turbine  as  a  function  of  head  and  flow  rate  can  be 
calculated  using  Eq.  (22).  Note  that  power  generated  is  a  linear  function  of  flow  rate  and 
head  available.  The  efficiency  for  the  turbine  given  is  approximately 


Pa,ht  =  opgQh. 


(22) 
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2.5.  Electrolyzer  modeling 


Electrolyzers  are  devices  that  use  electricity  to  split  water  into  hydrogen  and  oxygen.  The 
semi-empirical  model  of  an  alkaline  electrolyzer  was  reported  in  [63,64].  The  reactions  at 
the  anode  and  cathode  of  the  electrolyzer  cell  are  given  as  follows: 

Anode  :  H20  =  2H+  +  \02  +  2e~.  (23) 


Cathode  :  2H+  +  2e~  =  H2.  (24) 

The  semi-empirical  model  described  in  [63]  could  be  easily  applicable  to  PEM 
electrolyzers  with  a  change  in  the  empirical  constants.  The  voltage  required  for  the 
electrolyzer  is  obtained  from  the  summation  of  the  reversible  voltage  and  summation  of 
losses  due  to  various  resistances  in  the  cell.  The  losses  can  be  classified  as  activation  losses, 
which  quantify  for  slow  electro-kinetics,  ohmic  losses  in  the  cell  and  the  concentration 
losses,  which  occur  at  higher  current  densities  due  to  improper  mass  transfer  of  the  species: 


U ei  —  Urev-\-  rjact  +  rjohm  +  rj 


cone ' 


Reversible  voltage  as  a  function  of  temperature  and  pressure  can  be  given  from  the 
Nernst  equation  [64]  as  follows: 


1 

nF 


Uq  +  RT  In 


1/2 

Ph2Po2  , 


where  U0  is  the  standard  cell  potential  and  is  calculated  from  the  standard  Gibbs  free 
energy  using  Eq.  (27).  For  the  case  of  an  electrolyzer  reaction,  the  Gibbs  free  energy 
change  is  a  positive  quantity,  and  at  standard  conditions  (25  °C  and  1  bar)  it  is  equal  to 
237.13  kJ/mol: 


AG 

nF 


The  Gibbs  free  energy  is  a  function  of  temperature  and  pressure.  It  can  be  obtained  from 
a  change  in  the  enthalpy  of  formation  and  a  change  in  entropy  using  Eq.  (28).  The  change 
in  enthalpy  of  formation  and  entropy  as  functions  of  temperature  and  pressure  can  be  first 
obtained  for  individual  species  from  Eqs.  (31)  and  (32)  [65].  Then  the  other  parameters  can 
be  evaluated  for  the  electrolyzer  cell  using  Eqs.  (29)  and  (30): 


AG  =  AH  -  TAS, 


AH  =  AHHl  +  jAH0l  -  AHHl0, 
AS  =  A  Sh2  +  jASo2  —  ASh2o , 

T 

AH  =  AH°+  I  CpdT, 

7  E  ref 


(29) 

(30) 

(31) 


AS  =  AS0  + 


T 


Tref 


Cp 

T 


d T-R  In 


P 


\P  ref 


(32) 


The  specific  heat  at  constant  pressure  can  be  expressed  as  a  function  of  temperature 
using  Eq.  (33)  [66].  Constants  a,  /i,  y ,  3,  and  e  can  be  obtained  from  Table  2.  This  function 
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Table  2 

Constants  for  calculating  specific  heat  at  constant  pressure,  Eq.  (33)  [66] 


Gas 

a 

p  X  103 

y  x  106 

109 

£X  1012 

h2o 

4.070 

-1.108 

4.152 

-2.964 

0.807 

o2 

3.626 

-1.878 

7.055 

-6.764 

2.156 

h2 

3.057 

2.677 

-5.810 

5.521 

-1.812 

can  be  integrated  and  further  used  to  evaluate  changes  in  enthalpy  and  entropy  of 
individual  species  using  Eqs.  (34)  and  (35): 


C.,  =  R\a  +  pT  +  yT2  +  5T 3  +  e  Tl 


(33) 


PTZ  yT>  6T4  eT 


AH  =  AH°  +  R[  a T  + 


+ 


+ 


+ 


(34) 


AS  =  AS0  +  n(a  In  T  +  fiT  +  +  F-  + 

—  tflnf—  |.  (35) 

\Pref) 

The  reversible  potential  Urev  can  be  easily  determined  from  Eq.  (26)  as  a  function  of 
temperature  and  pressure.  Following  that,  the  electrolyzer  cell  potential  can  be  evaluated 
from  Eq.  (36).  Note  that  the  empirical  parameters  r,  s ,  and  t  can  be  obtained  by  curve¬ 
fitting  manufacturers’  data.  The  second  term  on  the  right-hand  side  quantifies  the  ohmic 
resistance: 

U el  =  Urev  +  —  I  ei  +  S  log  I  el  +  1^  ?  (36) 

where  the  empirical  parameters  r  and  t  can  be  defined  as  function  of  operating  temperature 
using  Eqs.  (37)  and  (38).  Note  that  the  temperature  is  in  Celsius  units  for  calculating  r 
and  t : 


r  =  ri+r2-^,  (37) 

d  ref 


t  —  t\  +  t2 


Tref 

T 


+  t3 


The  amount  of  hydrogen  produced  is  given  by 


nHl  =  Of 


N  c,ellel 

nF 


where  rjF  is  the  Faraday  efficiency  which  is  the  ratio  of  actual  hydrogen  produced  to 
theoretical  maximum  possible  production.  Similarly,  the  amount  of  water  required  for  a 
specific  amount  of  hydrogen  production  and  the  corresponding  amount  of  oxygen 
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produced  are  given  by 

n h2o  =  eh  2  =2  ho2-  (40) 

Performance  of  the  electrolyzer  is  dependent  on  the  operating  temperature.  The 
temperature  can  be  either  assumed  constant  or  can  be  calculated  by  using  a  thermal  model. 
A  lumped  thermal  capacitance  model  [63]  can  be  used  to  predict  the  temperature  of  the 
electrolyzer: 

d  T1 

^ t  ^  Qgen  Qloss  Qcoob  (41) 

where  Qgen  is  the  internal  heat  generation  of  the  electrolyzer,  Qioss  the  heat  lost  to  the 
ambient,  and  Qcooi  the  amount  of  cooling  required  to  keep  the  electrolyzer  at  a  rated 
condition.  These  are  given  as  follows: 

Qgen  =  ^ c,el(U ej  —  U  tn)I  eU  (42) 

where 


Utn  =  —=r9  (43) 

nF 

Q,oss=^iT-Ta).  (44) 

F t,el 

A  common  method  of  cooling  an  electrolyzer  is  to  flow  excess  water  through  the  unit  as 
compared  with  the  water  required  for  the  actual  hydrogen  consumption.  In  this  case,  it  can 
be  assumed  that  the  water  flow  related  to  the  cooling  requirement  is 

Qcooi  —  Ccw(Tcwj  —  T cw,o)  =  UAjlx  LMTD.  (45) 

It  is  worth  mentioning  that  there  is  a  paucity  of  publications  related  to  studies  of 
electrolyzers.  More  data  are  required  to  model  the  electrolyzer  properly  based  on  the 
physical  behavior.  It  is  imperative  to  provide  a  good  basis  for  the  semi -empirical 
parameters  used  in  the  model  for  calculating  /-  V  curve,  thermal  model,  etc. 

Using  the  above  approach,  the  effect  of  temperature  and  pressure  on  the  I-V  curves  of 
the  PEM  electrolyzer  is  plotted.  In  the  report  by  Cropley  [67],  experimental  data  were 
available  for  only  one  temperature  and  one  pressure.  Using  this  data,  the  empirical 
constants  were  obtained  for  the  best  fit  and  the  effect  for  different  temperature  and 
pressures  were  calculated.  See  Figs.  9  and  10. 

Some  recent  studies  have  focused  on  the  search  for  improved  and  cheaper  catalyst 
alternatives  for  hydrogen  and  oxygen  evolution.  The  performance  of  cobalt  and  nickel 
glyomines  is  found  to  be  inferior  as  compared  with  platinum  for  hydrogen  evolution  at 
cathode  [68].  Authors  reported  the  performance  of  Ir/Sn,  Ir/Ru-based  oxides  for  oxygen 
evolution  at  the  anode.  An  efficiency  of  94.4%  for  a  single  cell  PEM  electrolyzer  was 
achieved  with  Ir/Ru  at  anode  and  Pt/C  catalyst  at  cathode.  Ir/Sn  was  not  found  to  be 
effective  catalyst  for  oxygen  evolution  [69,70]. 
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2.6.  Fuel  cell  modeling 

A  fuel  cell  is  an  electrochemical  device  that  converts  hydrogen  and  oxygen  gas  into 
electrical  current  (DC).  A  review  of  fuel  cell  modeling  based  on  analytical,  semi-empirical, 
and  mechanistic  approaches  has  been  given  by  Cheddie  and  Munroe  [71].  Modeling  of  fuel 
cells  for  system  analysis  purposes  was  initially  proposed  by  Vanhanen  et  al.  [32], 
Chamberlin  [72],  Lee  [73],  and  Amphett  [74].  Further  improvement  in  the  fuel  cell  model  of 
Amphett  et  al.  was  reported  in  [75-79]. 

The  reactions  at  the  anode  and  cathode  are  given  as  follows: 

Anode  :  H2  =  2H+  +  2e~ .  (46) 

Cathode  :  2H+  +  ^02  +  2e~  =  FLO. 


(47) 
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Electrical  energy,  water,  and  heat  are  the  three  products  of  this  reaction.  The  total  cell 
potential  can  be  obtained  from 


U fc  Urev  T  Oact  "E  0 cone  "E  0 ohm’  (48) 

where  Urev  is  the  reversible  voltage  as  a  function  of  temperature  and  pressure  can  be 
obtained  from  the  Nernst  equation  [65].  The  first  term  gives  the  variation  of  standard 
potential  of  the  fuel  cell  with  respect  to  temperature  and  pressure: 


1 

nF 


Uojc  +  RT  In 


,PhJ>£ 


U0>fc  is  the  standard  cell  potential  of  the  cell.  It  can  be  calculated  as  a  function  of 
temperature  and  pressure  from  the  relations  developed  for  electrolyzers.  The  only 
difference  is  that  the  products  and  reactants  are  opposite  from  that  of  an  electrolyzer  and 
are  given  by  Eqs.  (50)  and  (51): 


AH  =  A HHl0  -  A HHi  -  \AH0l, 


AS  =  ASHl0  ~  ASh2  ~  \ASo2- 


Activation  polarization  r]act  is  the  extra  voltage  applied  to  counteract  the  slow  reaction 
rate  at  the  electrodes  of  fuel  cells.  The  reaction  rate  can  be  increased  either  by  increasing 
the  process  gas  flow  rate  and  thereby  increasing  the  process  gas  kinetic  energy;  or  by 
adding  extra  voltage  and  thereby  reducing  the  energy  output  of  the  fuel  cell.  According  to 
Faraday’s  law,  the  reaction  rate  is  proportional  to  the  current  density;  and  a  lower  reaction 
rate  implies  lower  current  densities.  Other  means  of  reducing  the  activation  potential 
besides  increasing  the  process  gas  flow  rates  include:  increasing  active  area  of  the  electrode, 
increasing  operating  temperature  and  increasing  the  utilization  of  catalyst  [73]. 

The  expression  for  the  activation  polarization  (anode  +  cathode)  was  derived  in  [78,79] 
from  Butler-Volmer  equation  as  follows: 


Oact  —  £i  +  In (cq2)  +  Q4T  ln(/), 

where  £1,  £2,  £3,  £4  are  given  as  follows: 


a  Ge; 
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Numerical  values  for  the  £l5  £2,  £3,  £4  parameters  were  provided  for  a  Ballard  fuel  cell 
[77,78]  as  follows: 


^  =  -  0.948(±0.004), 

£2  =  kDR  ^  +  k%  +  0.000197  In  G4/c)  +  4.3  x  105  In  (c*  2), 

^  =6.8  ±0.02  x  10“5, 

U  =  ~  1.97  ±0.05  x  10“4. 


These  values  could  be  changed  within  given  limits  for  modeling  PEM  fuel  cells  from 
other  manufacturers. 

Ohmic  polarization,  rjohm ,  is  the  extra  voltage  needed  to  overcome  the  electrical 
resistance  offered  by  the  cell  to  the  electron  and  proton  transfer  within  the  cell.  Note  that 
the  resistance  to  the  electron  transfer  could  be  due  to  electrodes,  graphite  collector  plates, 
and  other  fixtures  that  connect  the  membrane  electrode  assemblies  (MEA)  together.  This 
resistance  could  be  decreased  by  using  materials  that  are  good  electrical  conductors. 
Proton  conductivity  could  be  improved  by  reducing  the  thickness  of  the  membrane.  The 
ohmic  polarization  can  be  estimated  as  follows  [78]: 


d  ohm  i(R 


electronic 


+  Rpro,on )  =  -iR 


internal 


Here,  Relectromc  includes  resistance  due  to  contact  resistance.  For  pure  graphite  plates  the 
actual  Relectromc  should  not  be  significant  in  comparison  to  RProton  [79].  RProton  can  be 
obtained  as  follows: 


j^proton 


where  L  and  Afc  are  the  thickness  and  the  area  of  the  membrane,  respectively,  and  rm  is  the 
membrane  specific  resistivity  for  the  proton  transfer  from  the  membrane.  The  latter  is  a 
function  of  the  type  and  characteristics  of  the  membrane,  water  content  of  membrane, 
temperature,  and  current  density.  For  Nation  membranes  the  following  relationship  for  rm 
was  proposed  [78]: 


181.6 

[l  +0-03U)  +0-062  (3^3) 2  (4)  ] 

A0  +  Adr  age  —  0.634  —  3 

exp  (4.18  [T  j?03]) 

The  voltage  degradation  effect  on  the  performance  of  fuel  cells  was  studied  in  [79]. 
Aging  phenomena  are  introduced  in  the  model  through  the  A  terms  in  ohmic  losses  and  kceu 
in  activation  potential  losses. 

Concentration  polarization  rjconc  is  the  extra  voltage  required  to  overcome  the  resistance 
to  the  mass  transport  of  the  process  gases  at  the  reaction  sites  on  membrane  electrode 
assembly.  Mass  transport  includes  the  transport  of  gases  to  the  reaction  site  (flow  through 
channels  and  diffusion  in  the  electrodes)  and  removal  of  water  and  other  impurities  from 
the  reaction  site.  At  the  lower  current  densities,  these  losses  are  not  significant  as  compared 
with  higher  current  densities.  At  higher  current  densities,  process  gases  should  reach  the 
reaction  sites  at  higher  rates  and  similarly  the  water  and  other  impurities  should  be 
removed  more  quickly  from  the  reaction  sites.  Concentration  polarization  can  be  estimated 
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Current  Density  (mA/cm2) 

Fig.  12.  Effect  of  temperature  on  the  performance  of  a  fuel  cell. 


as  follows  [80]: 


rj 


cone 


RT 

nF 


lL 


Hence,  the  I-U  curve  of  the  PEM  fuel  cell  can  be  obtained  using  Eq.  (48)  considering  the 
losses  due  to  activation,  ohmic  and  concentration  polarizations.  Figs.  11  and  12  show  the 
performance  of  the  1.4  kW  fuel  cell  stack.  The  experimental  data  used  for  this  purpose 
were  reported  on  the  web  [81].  Figs.  11  and  12  show  the  effects  of  pressure  and  temperature 
on  the  performance  of  the  fuel  cell,  respectively. 

Similar  to  electrolyzers,  a  fuel  cell’s  performance  is  also  temperature  dependent.  In  a 
manner  similar  to  what  is  done  relative  to  an  electrolyzer  a  thermal  model  can  be  used  for 
predicting  the  temperature  of  the  fuel  cell.  Internal  heat  generation  for  this  case  would  be 
as  follows: 


Qgen  ~  Ncfc(Utn  Ufc)F 


(61) 
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In  addition,  a  fuel  cell’s  heat  loss  can  be  calculated  using  a  heat  balance  across 
the  inlet  and  outlet  streams  of  hydrogen  and  oxygen  flow.  In  the  case  of  a  hydrogen/air  fuel 
cell,  it  must  be  remembered  that  the  oxygen  side  balance  involves  21%  of  the  total 
atmospheric  air  flow.  Some  recent  studies  have  focused  on  accurately  predicting  the 
activation  and  concentration  polarization  [82-84].  The  evaluation  of  activation  and 
concentration  polarizations  at  the  anode  and  the  cathode  is  based  on  a  function  of  the 
concentration  and  partial  pressures  of  the  reactant  gases  (H2  and  02)  that  are  related  to 
each  other  by  Henry’s  law.  Once  the  exchange  current  density  is  obtained,  the 
approximations  of  Butler-Volmer  relationship  for  lower  and  higher  polarization  can  be 
used  to  obtain  the  activation  potential  for  the  anode  and  the  cathode  [83].  A  literature 
review  regarding  the  mechanisms  of  hydrogen  oxidation  reaction  on  different 
platinum  electrodes  and  their  effect  on  the  exchange  current  density  has  been  given  by 
Mann  et  al.  [84]. 

2. 7.  Storage  modeling 

Hydrogen  storage  is  one  of  the  important  aspects  of  the  path  towards  a  hydrogen 
economy.  As  was  noted  in  the  early  part  of  this  paper,  many  of  the  renewable  sources  have 
energy  outputs  that  are  quite  variable  in  time,  so  a  means  of  storing  the  results  of  capturing 
the  energy  is  critical.  The  following  are  the  commonly  used  hydrogen  storage  technologies 
[85-88]  currently  under  research: 

1.  Compressed  hydrogen. 

2.  Liquid  hydrogen. 

3.  Metal  hydrides. 

4.  Carbon-based  materials  (fullerenes,  carbon  nano  tubes,  activated  carbons). 

Discussions  for  all  these  forms  and  their  current  status  have  been  given  [85,87].  In  this 
section,  we  will  focus  on  modeling  aspects  of  these  technologies. 


2.7.7.  Compressed  hydrogen  storage 

Modeling  of  compressed  hydrogen  gas  storage  can  be  accomplished  by  assuming  with 
little  error  that  hydrogen  is  an  ideal  gas.  The  amount  of  work  required  to  compress 
hydrogen  to  a  required  pressure  can  be  obtained  by  considering  a  polytropic  process  for 
compressor.  The  actual  work  required  can  then  be  found  using  a  compressor  efficiency, 
information  usually  furnished  by  the  compressor  manufacturer,  using  Eq.  (62): 
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And  pressure  inside  the  storage  tank  can  be  obtained  from 

_  nstRT 
Pst  —  y 
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2. 7.2.  Liquefied  hydrogen 

In  general,  it  is  not  anticipated  that  this  approach  will  be  used  for  most  stand-alone 
hydrogen  systems.  Although  it  is  possible,  it  seems  that  the  nature  of  a  liquefaction 
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Fig.  13.  Shown  are  the  adsorption  curve  (upper  curve)  and  desorption  curve  (lower  curve)  for  a  metal  hydride 
hydrogen  storage  [11]. 

approach  may  be  too  complicated.  Of  course,  the  storage  energy  density  is  relatively  high 
using  a  liquid  compared  with  using  even  a  very  high-pressure  gas.  Good  sources  of 
information  on  this  general  topic  are  given  in  the  literature  [85,87,88]. 

2.7.3.  Metal  hydrides 

Many  metals  and  alloys  have  a  tendency  to  combine  with  hydrogen  to  form  hydrides 
[89].  When  hydrogen  is  stored  inside  a  hydride  (also  called  as  charging  or  hydriding),  the 
heat  is  released  during  this  filling  operation.  On  the  contrary,  heat  input  is  needed  while 
taking  hydrogen  out  from  the  hydride,  also  called  discharging  or  dehydriding.  This 
reversible  process  for  metal  M  can  be  expressed  as 

M  +  0.5x772  o  MHX  +  heat.  (64) 

Including  a  metal  hydride  with  the  renewable  hydrogen  system  needs  a  fairly  simple 
model  that  includes  the  rate  of  hydrogen  storage  (includes  hydriding  and  dehydriding),  the 
amount  of  heat  generated/required,  as  well  as  the  operating  temperature  and  pressure. 
Detailed  modeling  of  the  metal  hydrides  is  reviewed  in  [85].  A  database  of  metal  hydrides 
and  their  properties  is  available  on  the  web  [90].  Using  the  lumped  capacitance  method  the 
temperature  of  the  bed  as  a  function  of  time  can  be  calculated  [91]  as  well  as  the  heat 
generated  or  required  for  hydriding  or  dehydriding.  Based  of  bed  temperature,  the 
operating  pressure  can  be  found  using  Van  Hoffs  equation  [85]: 


Cmh  ^  =  Qmhit)  ~  keff[T(t)  -  T  a], 

(65) 

AH  AS 

(66) 

In  (Peg)  -  RT  R, 

where  Qmh  is  a  function  of  time  because  the  heat  would  be  released  during  the  absorption  and 
required  during  desorption.  A  basic  assumption  usually  made  in  the  modeling  of  hydride  storage 
is  that  the  cycling  efficiency  is  1 ,  which  means  that  there  is  no  change 
in  the  adsorption  and  desorption  properties  of  the  MH.  Fig.  13  shows  the  hydrogen  adsorption 
and  desorption  curve  at  5.3  °C  for  a  MH  bed  of  19  Nm  rated  hydrogen  storage  capacity  [11]. 
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Metal  hydride  beds  for  hydrogen  storage  integrated  with  RHGUS  have  been  reported  in 
the  literature  [11,30,34].  Experimental  data  for  the  metal  hydride  bed  integrated  with  an 
electrolyzer  and  a  fuel  cell  were  given  by  Vanhanen  et  al.  [92]. 

2.7.4.  Carbon-based  materials  for  hydrogen  storage 

Carbon-based  materials  are  capable  of  storing  hydrogen.  A  review  of  work  related  to 
hydrogen  storage  in  fullerenes,  carbon  nanotubes  and  activated  carbons  has  appeared  in 
the  literature  [86,93].  A  wide  range  of  hydrogen  uptakes  in  carbon  nanostructures  has  been 
reported  [93]. 

Fullerenes  are  synthesized  carbon  molecules  that  resemble  a  football  in  shape. 
Experiments  showing  adsorption  and  desorption  of  hydrogen  (up  to  6wt%)  have  been 
reported.  The  cyclic  behavior  of  these  materials  is  not  impressive  at  this  stage  (as  noted  in 
[86]). 

Two  types  of  carbon  nanotubes  have  been  explored:  single-wall  nanotubes  (SWNT) 
and  multi- wall  nano  tubes.  Hydrogen  can  be  stored  in  carbon  nano  tubes  by  chemisorption 
or  physisorption  mechanisms.  In  chemisorption  the  hydrogen  is  covalently  bonded 
with  host  materials.  The  desorption  temperatures  are  too  high  (>500K).  Physisorption 
occurs  due  to  Van  Der  Waaf  s  attractive  forces  between  hydrogen  and  the  host  material 
(carbon)  that  stabilizes  the  hydrogen  on  the  surface  of  carbon  nanotubes.  Desorption 
temperatures  of  hydrogen  for  SWNTs,  in  the  vicinity  of  room  temperature,  have  been 
reported  [94]. 

Activated  carbons  demonstrate  high  surface  areas  and  are  able  to  adsorb  hydrogen  in 
bulk  material  micro-pores  [86].  Experimental  data  reported  by  numerous  researchers  have 
been  tabulated  for  hydrogen  storage  (wt%)  in  carbon  nanotubes,  carbon  nanofibers — 
graphitic  planes  that  are  intercalated  at  a  fixed  distance  (0.355)  from  each  other — and 
activated  carbons  [93]. 

The  effective  integration  of  carbon  nanostructure  hydrogen  storage  with  RHGUS  or 
automobiles  will  require  the  development  of  this  material  with  a  capability  of  operating 
near  ambient  conditions  (25  <  T<  100  °C  and  P>  1  atm).  Effective  modeling  tools  must  be 
developed  to  predict  accurately  the  mechanisms  of  hydrogen  storage  in  carbon 
nanostructures.  The  cyclic  performance  of  carbon  nanostructures  needs  to  be  addressed 
along  with  the  degradation  of  performance  with  time. 

2.8.  Residential  application 

Energy  consumption  of  a  house  can  be  generally  categorized  in  two  forms,  electrical 
and  thermal  (as  a  fuel).  The  energy  demand  of  any  particular  house  is  a  direct  function  of 
the  type  of  materials  of  construction  and  their  ability  to  transfer  heat  from/to  the 
surrounding  atmosphere.  It  is  crucial  to  have  an  energy  efficient  design  for  future 
residential  houses  because  it  will  not  only  reduce  the  overall  energy  demand  of  the  house 
but  also  be  helpful  to  design  economically  feasible  RHGUS  configurations  for  this 
application. 

The  best  possible  scenario  for  designing  RHGS/RHGUS  for  this  application  is  to  use 
actual  energy  demand  data  for  the  house  on  an  annual  basis.  If  this  is  not  available,  then 
the  hourly  electrical  energy  demand  for  the  house  can  be  obtained  as  given  in  what  follows. 
Firstly,  the  cooling  load  and  heating  load  demand  can  be  estimated,  as  a  function  of  design 
conditions  and  type  of  materials  used,  for  the  particular  building  [48].  The  data  for  a 
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variety  of  structures  and  materials  and  their  R-values  (wall  structure,  roof,  basement, 
doors,  windows,  lights,  number  of  occupants,  etc.)  can  be  found  using  techniques 
summarized  by  ASHRAE  [95].  Secondly,  an  overall  heat  transfer  coefficient  for  the  house, 
based  on  floor  area,  can  be  estimated  for  the  design  conditions  of  a  particular  location. 
This  overall  heat  transfer  coefficient  can  be  further  used  to  calculate  the  hourly  electrical 
energy  demand  of  the  house  as  a  function  of  ambient  temperature.  As  an  example,  hourly 
energy  demand  for  a  standard  house  (1600  ft“)  located  in  Las  Vegas,  NV,  has  been 
reported  [48]. 

Once  the  hourly  electrical  energy  demand  is  known,  the  next  step  is  to  obtain  the 
thermal  energy  demand.  The  hydrogen-based  car  as  an  example  can  currently  go 
approximately  an  average  of  56.5  miles/kg  of  hydrogen  gas  [96].  Assuming  a  certain 
percent  of  home  refueling  (0-100,  where,  0  denotes  the  case  with  no  hydrogen  generated  at 
home  is  used  for  fueling  the  car  and  100  denotes  all  fuel  needs  of  car  are  satisfied  by 
hydrogen  generated  at  home)  can  give  the  hydrogen  required  for  driving  the  car.  Besides 
this,  the  gas  requirements  for  water  heating,  cooking,  and  drying  can  be  estimated  from 
data  given  by  EERE  [97].  Hourly  residential  gas  requirements  can  be  obtained  by  dividing 
the  total  amount  of  gas  required  by  the  number  of  hours  in  a  year. 

Recently,  the  potential  for  research  to  design  energy  efficient  buildings  has  been 
reported.  Numerous  aspects  of  energy  efficient  design  of  walls,  windows,  lighting,  etc.  have 
been  discussed  [98].  Another  study  reported  the  evaluation  of  evaporative-water  cooled 
air-conditioning  systems  as  compared  with  conventional  air-cooled  systems  for  a  house  in 
a  desert  climate  [99]. 

Based  on  these  electrical  and  gas  requirements  a  RHGS/RHGUS  system  can  be 
designed  and  its  performance  can  be  predicted  for  a  particular  house  for  any  location. 
Economic  issues  for  a  house  system  designed  according  to  the  information  given  above  are 
not  addressed  here.  However,  it  is  worth  mentioning  that  if  the  house  is  already  grid 
connected,  then  the  grid  can  be  easily  used  as  short-term  storage  as  compared  with  a 
battery.  Certainly,  this  would  not  be  a  case  for  an  off-grid  application. 

3.  Concluding  assessment  of  state  of  system  understanding 

From  this  review,  it  can  be  concluded  that  a  large  amount  of  effort  has  been  taken  in  the 
recent  past  to  model  the  various  components  of  RHGS/RHGUS,  either  based  on 
individual  components  or  collectively  for  the  overall  system.  Because  of  this  vast  literature 
on  simulation  and  experimental  efforts  the  question  that  might  be  asked  is  “Do  we  have  a 
good  understanding  of  the  methodology  to  model  the  individual  components  as  well  as  the 
overall  system  associated  with  a  RHGS/RHGUS?”  This  article  is  an  attempt  to  explore 
details  of  an  answer  to  this  question.  Although  simulation  information  presented  in  this 
paper  is  chosen  to  reflect  the  proper  physics  involved  with  the  major  components,  further 
research  is  needed  to  reduce/eliminate  generally  the  empirical  nature  of  the  models. 

The  input  energy  resources  needed  for  modeling  purposes  are  well-quantified  and  over  30 
years’  of  data  are  readily  available  for  many  locations  in  US.  PV modeling  discussed  here  is 
well  developed  and  can  be  applied  to  the  different  types  of  cells  currently  being 
manufactured.  Electrolyzer  models  are  less  well  developed.  For  example,  the  ones 
applicable  to  PEM  technology  are  currently  in  crude  form  and  need  attention  to  include 
more  of  the  nature  of  the  physics  involved.  Fuel  cell  modeling  has  received  several  years  of 
careful  attention,  and  this  area  is  fairly  mature.  Recently  published  papers  illustrate  that 
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general  performance  aspects  of  these  devices  can  be  predicted  quite  well.  The  physics  for 
the  description  of  hydrogen  storage  of  compressed  gas  are  well  understood.  Only  simple 
simulation  tools  have  been  put  forth  for  describing  metal  hydride  storage ;  more  data  are 
needed  on  large-scale  applications  to  model  accurately  this  method  of  containing  the  gas. 
Carbon-based  storage  technology  is  in  very  early  development  stages  so  its  ability  to  store 
effectively  on  large  scale  is  still  in  question.  On  the  other  hand,  residential  utilization  of 
hydrogen  can  be  readily  represented  in  simulations.  The  energy  requirements  for  a  house 
can  be  modeled  easily  using  methodology  that  is  well  developed.  However,  the  current 
designs  of  buildings  that,  in  general,  require  large  amounts  of  energy,  will  have  to  have 
large  RHGUS  installations,  and  these  will  be  very  expensive.  More  energy  efficient  design 
of  buildings  will  open  this  market  more  quickly. 

4.  Recommendations  for  future  work 

Based  on  the  conclusions,  some  of  the  key  aspects  from  a  modeling  standpoint  for  the 
RHGS/RHGUS  that  needs  to  be  addressed  by  the  scientific  community  are  listed  below: 

•  Application  of  RHGS/RHGUS  to  homes  is  limited  both  by  the  relative  immaturity  of 
the  technology,  as  well  as  the  relatively  large  systems  that  would  be  needed.  The 
development  of  standard  designs  of  energy  efficient  residential  houses  is  needed  to 
reduce  their  energy  demand.  This  will  assist  in  decreasing  the  sizes  of  energy  conversion 
equipment  substantially.  Appropriate  methods  to  include  the  energy  (electrical  and 
thermal)  efficient  equipment  and  their  effect  on  overall  energy  demands  needs  to  be 
quantified. 

•  Simulation  approaches  to  represent  the  performance  of  photovoltaic  cells  are  well 
developed.  Quantifying  the  degradation  of  performance  with  time  for  these  devices 
needs  to  be  addressed. 

•  Modeling  insights  for  electrolyzers  have  to  be  improved  significantly  to  show  clearly  the 
activation,  ohmic  and  concentration  polarization  losses  inside  the  cell.  Since  there  is  not 
enough  data  available  for  PEM  electrolyzers,  these  aspects  certainly  should  be 
addressed.  Models  need  to  be  developed  so  that  actual  mechanisms  of  electrolysis  can  be 
reflected  in  it  instead  of  theoretical  estimates. 

•  CFD  modeling  has  been  widely  used  to  simulate  the  multi-physical  phenomena 
associated  with  fuel  cells.  Not  nearly  as  much  attention  has  been  given  to  the  PEM 
electrolyzers,  and  the  understanding  of  this  technology  could  gain  a  great  deal  from 
these  types  of  efforts. 

•  Work  related  to  the  hydrogen  storage  technologies  has  been  reported  in  [85].  The 
physics  of  some  aspects  of  this  important  component  of  renewable  systems  are  well 
understood,  but  there  are  several  newer  approaches  that  may  hold  great  promise  that 
need  to  be  investigated  further. 

•  The  degradation  of  performance  validated  with  actual  experimental  data  needs  to  be 
addressed  for  almost  all  the  components  of  RHGUS. 

•  Since  the  technology,  manufacturing,  and  marketing  of  many  of  the  RHGUS 
components  are  in  early  stages  of  development,  costs  of  all  the  components  of  the 
RHGUS  are  very  uncertain.  This  topic  should  be  addressed  now  with  prices  estimated 
for  the  situation  when  this  technology  becomes  more  mature,  and  then  revisited 
periodically  as  the  market  continues  to  grow. 
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